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ABSTRACT 


This thesis investigates intermodulation distortion generated by analog-to-digital 
converters (ADCs) in a channelized digital ESM receiver when processing multiple signals 
simultaneously. Spurious free dynamic range (SFDR) associated with this distortion is 
discussed. Two methods for increasing spurious free dynamic range are evaluated. First, by 
adding a small amount of Gaussian noise to the input of the receiver, the inteimodulation 
distortion is found to be reduced significantly. Second, by using a narrow bandwidth sub- 
Nyquist sampling rate with high dynamic range ADCs it is possible to increase the spurious 
free dynamic range of the digital receiver. The first method is a simple approach but the 
ability to increase the SFDR is limited. The second method is more effective but requires 
greater computation and complex receiver design. 
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L INTRODUCTION 


A. CHANNELIZED DIGITAL ESM RECEIVER AND ANALOG-TO-DIGITAL 

CONVERTERS (ADC) 

Modern electronic support measures (ESM), or electronic support (ES), receivers 

I 

require a wide instantaneous bandwidth and high dynamic range in order to receive and 
process multiple signals at the same time. Modem signals such as frequency agility radar may 
have a 500 MHz frequency hopping bandwidth or greater. It is impossible for a narrow band 
superheterodyne ESM receiver to tune to each frequency component of the frequency agility 
radar. However a channelized digital ESM receiver with a 500 MHz or higher instantaneous 
bandwidth using a fixed local oscillator frequency, can process all frequency components 
within its bandwidth by taking samples of in phase and quadrature sinusoid components from 
the output of an analog-to-digital converter (ADC) and utilizing the Fast Fourier Transform 
(FFT) to perform frequency domain analysis. This is accomplished without rapidly tuning the 
local oscillator frequency as occurs in a conventional narrow band superheterodyne receiver. 

Spectral analysis techniques are the key to advanced receiver designs and provide an 
essential parameter for emitter sorting and also for the efficient allocation of ECM resources. 
Real-time digital spectmm analysis is possible using special hardware configured to 
accomplish the fast Fourier transform (FFT) algorithm. High speed multipliers and 
accumulators are currently available in large-scale integration (LSI) form, which can be linked 
together in pipeline structures to perform a complex FFT. A critical component which limits 
the use of a high speed digital processor is the analog-to-digital converter (ADC), which 
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translates the analog waveform at the input of the processor into a digital word that is 
manipulated in the processor. [Ref 1] 


Due to the rapid development of analog-to-digital converters (ADCs), building of 
wideband digital ESM receivers becomes a reality. Digital ESM receivers rely on high speed 
ADCs and digital signal processing techniques in order to encode the received signals into 
Pulse Description Words (PDW) in real time and store them into memory for future analysis 
and ECM functions. Figure 1 shows a survey on the capability of current ADCs. The 
information on the ADCs was obtained from [Ref 2] through [Ref 10]. 
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Figure L Performance of Current ADCs 


A channelized digital ESM receiver can be represented in functional blocks as shown 
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in Figure 2. The output from the ADC is digital. These data are in the time domain and the 
information is available as spectral lines or spectrum density. However, these outputs do not 
completely satisfy ESM requirements. The spectral lines must be converted into carrier 
frequencies of the input signals. In order to emphasize this process, a parameter encoder is 
identified separately from the spectrum estimator. The parameter encoder converts the 
frequency information into the desired Pulse Description Words (PDW) [Ref 10]. 



Digital EW Receiver 


Digital 

Processor 


Figure 2. Block Diagram of a Typical Channelized Digital ESM Receiver 


Since the channelized digital ESM receiver allows multiple radar signals to be down 
converted into IF components within the passband at the same time, intermodulation products 
appear at the ADC output and the receiver reports more signals than it has received. We can 
generally increase the threshold setting of the receiver to eliminate intermodulation products, 
but these will decrease the dynamic range of the receiver and lose weak signals. Sometimes 
the weak signals are the most important signals because they can come from weapon control 
radars or airborne fire control radars that are transmitted at low power levels intentionally to 
reduce the probability of being intercepted. The intermodulation products created by the 
ADC do not always behave as they do for linear devices where reduced input levels result in 
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predictable reductions in intermodulation products. We need to find other ways to reduce the 
intermodulation products and increase dynamic range to an acceptable level which is the 
objective of this thesis. 

B. THE EW ENVmONMENT AND PULSE DENSITY 

According to the research done by Peot [Ref 11], the worst case radio frequency 
(RF) electronic combat environment is a high-altitude aircraft conducting a mission over 
central Europe. Table 1 shows the total emitter and pulse density in the year 2000, which 
includes the signals from airborne intercept radars, ground-based SAM, ground-based AAA 
and other miscellaneous emitters. These figures are based on rough assumption in the year 
2000 which may differ from the actual environment. The assumptions for the scenario are as 
follows: 

• No more than 25% of the aircrafts are in the air at once for either side and 50% 
of the airborne intercepter radars are on all of the time. 

• The average pulse repetition frequency (PRF) for ground-based or miscellaneous 
emitters is 200 Hz. 

• The PRF for airborne fire control radar is 15 KHz. 

• The receiver is assumed operating in noiseless environment. 

• The terrain is flat. 

• About 50% of the AAA system will be radar-guided. 

From Table 1 we can see that there are several million samples per second for each 
type of emitter from each side. It is the number of pulses impinging on the aircraft. The 
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actual number of pulses received by the aircraft is scaled by the receiver sensitivity, terrain 
masking, receiver antenna coverage, look through percentage and other factors. As stated 
by Peot, using parallel architecture and high speed special purpose processing the EW system 
designer can achieve the processing densities required for the EW signal processors of the 
year 2000. 

The probability of pulse overlaps will be high in the above scenario. A good 
channelized digital ESM receiver design should be able to reduce the intermodulation 
distortion which may occur in such environment. 


Table 1. Total Emitter and Pulse Density Figures 


Emitter 

No, Of 
Emitter 

% in 
LOS 

% 

Active 

Emitters 

Visible 

Mainlobe 

Sidelobe 

Visibility(%) 

Average 

PR 

Pulse 

Density 
Contribution 
(Pulse per second) 

NATO 

A/A Radar 

1250 

100 

50 

625 

50 

15000 

4.7x 10« 

Warsaw Pact 
A/A Radar 

2000 

100 

50 

1000 

50 

15000 

7.5 X 10' 

NATO 

AAA 

4250 

50 

100 

2125 

50 

2000 

2. lx 10® 

Warsaw Pact 
AAA 

2650 

50 

100 

1325 

50 

2000 

1.3 X 10® 

NATO 

SAM 

4400 

50 

100 

2200 

50 

2000 

2.2 X 10® 

Warsaw Pact 
SAM 

19000 

50 

100 

9500 

50 

2000 

9.5 X 10® 

Misc 

1000 



1000 


2000 


Total 

34550 

— 

— 

18000 

— 

1 

1 

1 

29x 10® 
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n. FREQUENCY DOMAm MEASUREMENT OF INTERMODULATION 
DISTORTION AND SPURIOUS FREE DYNAMIC RANGE 

Channelized digital ESM receivers employ baseband and direct IF digitization which 
place stringent performance requirements on analog-to-digital converters (ADCs). Receiver 
performance is limited by the dynamic parameters of the ADC. We can derive the dynamic 
parameters using the Fast Fourier Transform (FFT). The MATLAB software provides unique 
FFT and FFTSHEFT functions which are convenient to use in digital receiver simulation. All 
the simulations and calculations in this thesis were done in MATLAB. 

The input signal is set to be a IV peak amplitude sine wave. The ADC simulation 
converts the IV peak amplitude sine wave signal to its maximum response which means the 
maximum output of the ADC is also set to 1 V peak. Any analog input signal is rounded (or 
truncated) to its nearest digital code, the rounding (or truncating) error represents the 
quantization error in a real ADC. The analog input range is divided into 2" digital codes so 
that each digital output code represents a corresponding analog input level within the 
quantization range. The FFT operation collects the digital codes and generates data for 
frequency domain evaluation. 

The key frequency domain ADC specifications concerning this thesis are signal-to- 
noise-ratio (SNR), total harmonic distortion (THD), dynamic range (DR), intermodulation 
distortion (IMD), signal-to-noise plus distortion (SINAD), and spurious free dynamic range 
(SFDR) which can be obtained from the FFT operation. 

The most important dynamic specification of a converter is the signal-to-noise ratio 
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(SNR). The SNR is related to the quantization noise power of the ADC, which determines 
the sensitivity level of the receiver. The formula for calculating the maximum SNR for an 
ideal ADC is 6.02n+1.76 dB where n is the number of quantization bits [Ref 10]. SNR 
measures signal power relative to noise power which is appropriate to single input signal 
analysis. We can calculate the actual noise power from the FFT output. SNR as applied to 
a single signal would be calculated as 

SNR=l0log(^^) dB (1) 

where is the power of the input signal after FFT processing and is the sum of all 
quantization noise after FFT processing, or, equivalently, by removing the carrier and 
harmonic components from normalized FFT spectrum (set the ideal full-scale spectrum level 
to 1) and calculate the RMS values of the remaining points. Since the SNR is in positive 
decibels, negative sign is needed. We can simplify the Equation 1 as: 

SNR^-mc,g(lP^J dB (2) 

For example, if there are 512 samples in the normalized FFT spectrum, then we can 
remove the carrier sample and sum up the remaining 511 quantization noise samples to obtain 
the SNR without considering the carrier amplitude since it has been normalized to 1. The rest 
of the dynamic specifications may follow this principle. 

THD measures the total power at the harmonic distortion frequencies Pharmomcs by 
removing the carrier and noise components from the normalized FFT spectrum and 
calculating the RMS values of the remaining points. THD is given by the following form: 
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rao=ioiog(EP ,dB 


(3) 


IMD measures the total power at the intermodulation distortion frequencies by 
removing the carrier signals, harmonics, and noise from the normalized FFT spectrum and 
calculating the RMS values of the remaining points. IMD is given by the following form: 

/A/Z)=101og(£P^„J dB (4) 

SINAD measures the total power at the harmonic distortion frequencies and noise by 
removing the carrier signals from the normalized FFT spectrum and calculating the RMS 
values of the remaining points. This is more appropriate when analyzing multiple input 
signals. SINAD is given by 

intermods ) dB (5) 

There are three different types of dynamic range we need to consider for evaluating 
the performance of an ADC; 

• Single signal dynamic range (DR): 

DR measures the power ratio of the strongest signal that the receiver can 
properly detect without generating spurious responses to the signal at the 
receiver’s sensitivity level. This is equal to the signal-to-noise ratio (SNR) of 
the system measured over a bandwidth equal to half the sampling frequency. 

• Spurious free dynamic range (SFDR): 

SFDR relates to the third-order intermodulation or harmonic distortions which 
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measures the power ratio of the strongest signal that the receiver can properly 
encode without generating detectable third-order intermodulation or harmonic 
distortions to the power at the noise level. It is equal to the ratio between the 
maximum signal component and the largest distortion component that can be 
obtained. 

• Instantaneous dynamic range (IDR); 

IDR measures the power ratio of the maximum and minimum simultaneous 
received pulses that can be properly encoded by the receiver. 

The three different types of dynamic ranges can be calculated from the normalized 
results after FFT processing by using the following equations; 

DR^Sm ( 6 ) 

SFDR = -1 Olog (Peak Harmonic or Intermod Power) (7) 

IDR =-1 Olog (Peak Weakest Signal Power) (8) 


A. SIMULATION OF A CHANNELIZED DIGITAL ESM RECEIVER 

Figure 3 shows the basic simulation structure of a channelized wideband digital ESM 
receiver using the two-channel approach. The passband (Hgp(f)) of the receiver is designed 
from 1000 MHz to 1500 MHz which can detect the lower-half of L-band search radar signals. 
The input signal passes through a bandpass filter and a 90 degree hybrid divides the input 
signal into in-phase (I channel) and quadrature-phase (Q channel) channels. The local 
oscillator generates 1250 MHz sine wave signal that feeds directly to the I channel and the 
Q channel mixers, The mixers down convert the input signals into two baseband video 


10 



signals. The I and Q signals feed two separate lowpass filters. The cutoff frequency of the 
lowpass filters is 250 MHz. The bandwidth of the receiver equals the sampling rate 

of the I-Q channel ADC. The ADCs digitize the IF signals into binary codes with maximum 
response to the input signal. Since the I-Q outputs are both digitized, the information content 
of the input signal is essentially doubled. 



Figure 3. Block Diagram of a Two-Channel Digital ESM Receiver 


For a N point FFT with sampling rate equal to half of the Nyquist rate/,=/^2, the 
frequency resolution AF for each FFT bin would be: 

A 

AF=^ (9) 

N 

The bandwidth (BW) of a two-channel digital ESM receiver equals the sampling rate f;. 

BW-l ( 10 ) 

The N point FFT functions as N bandpass filters, each FFT bin covers a frequency band as 
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shown on Figure 4. 



B. SINGLE SIGNAL ANALYSIS 

The simulation is done by using a 1 V peak amplitude sine wave signal sin(2;5^^) as 
an input signal to the receiver. The in-phase local oscillator signal is lsm.{27rfj) which 
directly feeds to the I channel mixer, while the Q channel mixer signal is 2sin(27z/;^+90°) or 
IcosilTrfJ). In order to simplify the simulation, we only consider the frequency components 
below 250 MHz and neglect the effects of the lowpass filters. The outputs from the I and Q 
channel become: 

7(0=sin(27iy;,0* 2sin(27r/p0 

=[cos(27t/-2Tt/o)t-cos(2iTy;.+27t/o)t] (11) 

IjJt)=cos{2Tzf-2Tif^i 

Q(t) =sin(27ty;0* 2cos(27r/Q0 

=[sin(27t};-2iT/o)?-sin(27ty;.+2Tt/(,)/] (12) 

Q^{t)=sm{2nf-2nf,)t 

The high frequency components sm{27f,+2 nf„)t and cos(2 7tf,+2 Trfjt are filtered out 
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by the lowpass filters resulting in the last lines of Equation (11) and (12). The lip(t) and 
Qi/t) are sent to two separate ADC which after digitization become I(n) and Q(n). The l(n) 
and Q(n) are combined as I(n)+jQ(n) for the FFT process. 

Figure 5 shows the output samples from the I and Q channels after an 8 bit ADC. The 
input 1210 MHz sample signal is down converted to -20 MHz IF signal. The number of 
quantization levels of the ADC isM=2"'‘=2*'*=128 (8 bit ADC is 7 bits plus sign) for the 
positive and negative cycles of the IF signals. After executing the FFTSHIFT function in 
MATLAB, the frequency information will be shifted back to 230 MHz. 



Time (ns) 


Figure 5 . ADC Output of a Digital Receiver with One Input Signal 

Figure 6(a) shows a normalized EF spectrum of a digital ESM receiver using an 8 bit 
ADC and 512 sample FFT with one input signal at 230 MHz. The sampling rate for the ADC 
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is set to the Nyquist rate of 500 MHz. The data shown from sample 1 to sample 512 is for 
500 MHz IF bandwidth. Since the amplitude of the carrier is a function of the converter’s 
resolution in bits, normalizing the amplitude of the main frequency component to set the ideal 
full scale input level to 0 dB simplifies the evaluation and dBc means dB below carrier 
amplitude. We can use the normalized FFT results to find the peak noise level and determine 
the dynamic range of the receiver for a single input signal. Since there is only one input 
signal, there is no intermodulation distortion, but only the spurious response and quantization 
noise. Figure 6(b) shows the comparison of theoretical SNR and simulation results as a 
function of the number of ADC bits for one input signal. The simulation program for single 
signal analysis is listed in Appendix A. 



ADC Bits 


Figure 6 . IF Spectrum of a Digital Receiver with One Input Signal 
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Table 2 shows a comparison of the theoretical one input signal SNR and a computer 
simulation with results for different numbers of ADC bits for 230 MHz input signals. The 
simulation results for 230 MHz input signal are within 3dB from the ideal SNR values and 
increase at 6dB per bit. 


Table 2. Theoretical SNR and Simulation SNR vs. Number of ADC Bits 


ADC Bits 

Ideal SNR (dB) 

Simulation SNR (dB) 

1 

7.78 

6.34 

2 

13.80 

12.02 

3 

19.82 

17.76 

4 

25.84 

23.74 

5 

31.36 

29.69 

6 

37.88 

35.59 

7 

43.90 

41.74 

8 

49.92 

47.76 

1 

9 1 

55.94 

53.41 

10 

61.96 

60.37 

11 

67.98 

65.85 

12 

74.00 

72.68 


C. TWO SIGNALS ANALYSIS 

Figure 7(a) shows the normalized output spectrum of a digital receiver with two 1 V 
peak amplitude sine wave signals/;=1230 MHz and^^=1270 MHz as the test signals. We can 
clearly see from the plot that there are ten equally separated inband intermodulation products 
generated by the two strong input signals which are the third { 2 f,-f 2 and 2/j-/;), fifth 
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and 3fr2f^, seventh { 4 f,- 2 f 2 and •?//), ninth {, 5 f,- 4 f 2 and 3 f 2 - 4 fJ, and eleventh and 

6 f 2 - 5 fJ order intermodulations whose amplitudes are higher than the peak noise level. 

Figure 7(b) shows the comparison of ideal SNR and simulation results for different 
numbers of ADC bits. 



ADC Bits 


Figure 7. IF Spectrum of a Digital Receiver with Two Input Signals 


Because of the intermodulation products, the SINAD obtained in this simulation is 
significantly worse than the theoretical SNR and the simulation result for a single signal SNR 
in the last section. The simulation program for two signals analysis is listed in Appendix B. 

The intermodulation generated by the ADC results from the nonlinear quantization 
effects. Quantization introduces not only noise and distortion of the input signals, but also 
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generates intermodulation products from the two strong input signals. When the frequency 
spacing between f, and £ is smaller than the bandwidth of the receiver, the third order 
intermodulation becomes the principal source of spurious in band signals since it has the 
highest amplitude of all intermodulations. Figure 8 shows the quantized output samples of 
the I and Q charmel ADCs for two input signals. 



Time (ns) 

Figure 8 . ADC Output of a Digital Receiver with Two Input Signals 


The dynamic range in this case refers to the two-signal or multiple signals analysis and 
is called the spurious free dynamic range (SFDR). SFDR is the most important specification 
in digital receiver design which determines the lowest distinguishable level of the input signals. 
As opposed to the sum of signum function used by Blachman [Ref 14], the quantization 
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staircase output is represented by the sum of the input plus a sawtooth. Using an Fourier 
series for expanding the sawtooth, series expressions are obtained for the output quantized 
carriers and the resulting intermodulation products of any order. Blachman obtained an 
expression for the amplitude of the third-order and higher order intermodulation due to 
quantization. Because of the diflficult mathematical integrations, only some special cases have 
been analyzed and further studies are required. Due to the complexity of Blachman approach, 
the procedure is not be followed in this thesis. 

Table 3 shows the simulation results of SFDR, SINAD, and MD for two input signals 
and verses the number of bits in the ADC. As expected, we can see that the intermodulation 
distortion decreases when the number of ADC bits increases. 

Figure 9 shows two strong signals (210 MHz and 270 MHz) accompanied by a weak 
si g nal (230 MHz) whose amplitude appears between the third-order intermodulation and the 
peak noise level. The receiver does not recognize this weak signal as an input signal and 
tends to ignore it. Chapters III and IV will discuss methods of reducing intermodulation 
distortion and restoration of the weak signals. 
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D. MULTIPLE SIGNALS ANALYSIS 


When multiple signals are applied to the input of an ADC, the output spectrum 
becomes unpredictable. Some intermodulation products tend to cancel with quantization 
noise resulting in amplitude reduction. The relationship between carriers and intermodulation 
products become uncertain when multiple signals are applied simultaneously to the receiver. 
The larger the number of input signals the worst the situation. In a real ESM environment, 
the reports from the receiver will be unreliable. Figure 10 shows the output spectrum of the 
ADC with fifteen sample signals and 8 bits quantization. 



Figure 10. IF Spectrum of a Digital ESM Receiver with Fifteen Input Signals 
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When multiple signals are present, the number of intermodulation products becomes 
unpredictable and much smaller in amplitude when compared to Figure 10. Figure 10 clearly 
shows the frequency spacings between carriers and each intermodulation product as well as 
the amplitude of each intermodulation product. Table 4 lists the SFDR and SINAD of an 8 
bit ADC for different numbers of input signals. We can see from the table that the SFDR and 
SINAD decrease as the number of input signals increase. As the number of input signals 
reached 15, the SFDR is almost equal to the SINAD. 


Table 4. SFDR and SINAD of an 8 bits ADC 


# of Signals 

SFDR (dB) 

SINAD (dB) 

2 

51.50 

40.10 

3 

43.67 

41.43 

4 

45.07 

40.82 

5 

45.55 

39.99 

6 

48.04 

39.14 

7 

45.02 

38.50 

8 

41.81 

37.57 

9 

42.24 

37.42 

10 

39.37 

36.56 

11 

39.00 

36.27 

12 

37.61 

35.58 

13 

38.57 

35.53 

14 

37.34 

35.00 

15 

34.72 

34.74 
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m. REDUCTION IN INTERMODULATION DISTORTION DUE TO ADDITIVE 

NOISE 

When designing an analog receiver, the noise level is kept as low as possible in order 
to maintain a higher receiver sensitivity. In a digital ESM receiver, which uses an ADC as a 
crucial component, it is not necessary to follow this rule. Noise sometimes has a positive 
effect in an ADC as proposed by Tsui [Ref 10] and Blachman [Ref 15]. Noise may reduce 
the spurs generated by quantization. In fact, the effect of the noise is to smooth the nonlinear 
characteristic of the ADC. In general, the signal builds up by a factor of N for an N-point 
FFT with respect to noise. The FFT process has a gain of lOlog (N) dB for the coherent 
input signals and 0 dB for the noncoherent noise so adding noise to the digital receiver does 
not affect the amplitude of the carrier. Figure 11 shows the block diagram of a digital ESM 
receiver with a noise source added after the bandpass filter. 



Figure 11. Block Diagram of a Digital ESM Receiver with a Noise Source 

The digitization coherence will be reduced by the noise power as well as the 
intermodulation products. The intermodulation products may be reduced or even disappear 
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by controlling the additive noise level properly, but excess input noise power also reduces the 

SINAD level and reduces the SFDR of the digital ESM receiver. 

The effect of adding noise to the input upon the distortion and intermodulation 

depends only on the probability distribution of the noise and is independent of its spectrum 

[Ref 15]. The presence of a small amount of Gaussian wideband noise which is readily 

eliminated by a low pass filter, should cause the intermodulation products due to quantization 

to be undetectable. 

A. ADDITIVE NOISE LEVEL AND INTERMODULATION DISTORTION 

Figure 12 shows the comparison of the IF spectrums of the two input signals with and 
without additive noise. The two input signals are/;=1230 MHz andy^=1270 MHz and the 
ADCs use 8 bits. In Figure 12(a), we can observe the intermodulations between// and/. 
In Figure 12(b), we can see the amplitudes of the intermodulation products decrease resulting 
from -40 dBc additive noise. When adding noise to the input, the SFDR increases but also 
causes SINAD to decrease. Its a trade off between SFDR and SINAD. 

There is no mathematical approach to obtain the optimum additive noise level for 
reducing the intermodulation products. In computer simulation, we can calculate the 
intermodulation distortion for different additive noise level and find the desired setting, but 
in practice we need to manually adjust the noise levels during intercepts. 
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Figure 12. The Out Spectrum (a) with or (b) without Additive Noise 


Figure 13 shows the relationship between additive noise level and intermodulation 
distortion of an 8 bit ADC for two input signals. The additive noise power starts from -140 
dBc to 0 dBc for testing of two sample signals. When noise is added, the intermodulation 
distortion appears smaller until the noise power reaches -40 dBc, then all the intermodulation 
products become noise. The highest additive noise power for this simulation should be -40 
dBc; levels higher than -40 dBc will have a negative effect on the receiver. 

Figure 14 shows the relationship between the additive noise level and SINAD of an 
8 bit ADC for two input signals. When noise is added, SINAD starts to decrease until the 
noise power is greater than -40 dBc, then the SINAD starts to drop drastically. 
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B. ADDITIVE NOISE LEVEL AND SFDR 


Figure 15 shows the relationship between additive noise level and SFDR of an 8 bit 
ADC for two input signals. When noise is added, the SFDR starts to increase until the noise 
power reaches -40 dBc, then the SFDR declines drastically. In comparison with the 
relationship between additive noise and intermodulation distortion and SINAD, the additive 
noise is most beneficial to the SFDR. 



Figure 15. Additive Noise level verses SFDR 


In order to test the SFDR after adding noise to the mixers, we can add a third input 
signal to the receiver to see if it is possible to restore the third signal from the intermodulation 
distortion. A third input signal/j=1240 MHz is used to perform the test. The ADC is still 
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8 bits. The amplitude of the third input signal is set to -50 dBc for testing. We can see from 
Figure 16(a) that when a third input signal is added to the spectrum, it generates more 
spurious signals rather than the intermods. 

By adding -40 dBc noise to the input, we can examine the amplitude of the third input 
signal which becomes higher than the third-order intermodulation products that are generated 
by/; and/ and all of the spurious signals disappear. Figure 16(b) shows that the amplitude 
of the third input signal now becomes greater than the third-order intermodulation products 
by the increased SFDR of the receiver. Again, there is no mathematical formula to specify 
the relationship between additive noise level and SFDR, but the noise does have a positive 
effect which increases the SFDR of the receiver by a small amount. 




Frequency (MHz) 

Figure 16 . SFDR Test with a Third Input Signal 
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rV. SUB-NYQUIST RATE APPROACH 


In the previous chapter, we have seen how additive noise has the effect of reducing 
intermodulation distortion and increasing the SFDR, but the SFDR is still limited by the 
number of ADC bits unless we narrow the IF bandwidth in order to use an ADC with a larger 
number of bits to achieve higher dynamic range for the digital ESM receiver design. 
Narrowing the IF bandwidth means more receivers are needed in the system which makes the 
system more expensive and difficult to maintain. 

An effective way to obtain a wide bandwidth with low sampling rate and larger 
number of ADC bits is to use sub-Nyquist rate sampling. The sub-Nyquist sampling scheme 
is very similar to the instantaneous frequency measurement (IFM) receiver. They both 
correlate the delayed and undelayed channel to obtain the frequency information. 

Figure 17 shows the basic structure of a digital ESM receiver using the sub-Nyquist 
sampling approach. The digitized IF outputs are processed by the FFT operation. The FFT 
operation generates real and imaginary parts for the delayed and undelayed channels in the 
frequency domain. 

Let XJk) and X^Jk) represent the real and imaginary parts of the undelayed channel, 
?tnAXf/k) be the delayed channel, and x is the delay time. The amplitude information 
can be calculated from the undelayed channel FFT output as: 

(13) 

and the frequency information can be calculated from the phase difference between the 
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delayed and undelayed channels. 



Figure 17. Basic Sub-Nyquist Sampling Structure of a Digital Receiver 

Let 6„ be the phase of the undelayed channel and 6^ be the phase of the delayed 
channel. The frequency information can be calculated as follows: 


0, =tan T-] 

“ "xjkY 

(14) 

0.=tan-n ] 

(15) 

A0=0^-0„=2tct 

(16) 

0^-0 

d u 

Inx 

(17) 
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The delay time x determines the highest frequency as well as the bandwidth of the 
receiver (not the sampling rate/,). When 0j-6^ equals to 2 k, the result in Equation (17) will 
be the highest frequency component from the FFT operation which represents the bandwidth 
of the receiver and has the following form: 

AD r 1 

=- ( 18 ) 

Figure 18 shows the results of the sub-Nyquist rate simulation. The two input signals 
are 1230 MHz and 1270 MHz, the IF bandwidth is still 500 MHz, but the sampling rate is 
reduced to 100 MHz (1/10 Nyquist rate) to accommodate the sub-Nyquist rate scheme. The 
delay time setting for the delayed channel is 0. Ins. The maximum unambiguous bandwidth 
for the channelized digital ESM receiver can be derived from Equation (18) which is 
AB=l/0. lns=10 GHz. The sampling rate has been reduced to 100 MHz, thus, we can select 
a 12 bit ADC with 100 MHz bandwidth instead of using an 8 bit ADC with IGHz bandwidth. 
The desired 500 MHz IF bandwidth will fit entirely into the 100 MHz bandwidth of the 12 
bit ADC without ambiguity. An 8 bit IGHz bandwidth ADC is also more expensive and has 
lower dynamic range than a 12 bit 100 MHz bandwidth ADC. 

The input signals have been down converted to 230 MHz and 270 MHz by the mixers. 
From Figure 19, the IF signals are measured at 30 MHz and 70 MHz in the 100 MHz wide 
IF spectrum. Since the corresponding unambiguous frequency band of the channelized digital 
ESM receiver is 10 GHz, we do not know the exact frequencies of the two input signals 
because it can be 30 MHz and 70 MHz, 130 MHz and 170 MHz, 230 MHz and 270 MHz, 
330 MHz and 370 MHz, 430 MHz and 470 MHz,..., or 9930 MHz and 9970 MHz. By using 
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the phase calculation (Equation 17), we can obtain the true frequency information of the input 
signals correctly without ambiguity. A simulation program for the sub-Nyquist sampling 
scheme is listed in Appendix C. 

The intermodulation distortion reduces when the input signals are the same as the 
previous simulation that used the Nyquist sampling rate, but intermodulation distortion does 
occur for some input signals. By adding noise to the input of this scheme, the intermodulation 
distortion can be reduced to an insignificant level. 


(a) Output Spectrum of a 12 Bits ADC Sampling at fs=100MHz Sub-Nyquist Rate 
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(b) Output Spectrum of an 8 Bits ADC Samplimg atfs=500MHz Nyquist Rate 



Frequency (MHz) 


Figure IS.The Comparison of Sub-Nyquist and Nyquist Sampling Spectrum 


Figure 19 shows the results when a 1240 MHz input signal is added to the receivers 
using a different number of bits for the ADCs. The amplitude of the third input signal is set 
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to -70.56 dB below the strongest amplitude of the carrier. Figure 19(a) shows a receiver 
using a 12 bit ADC with sub-Nyquist sampling rate; it is possible for the receiver to achieve 
a SFDR greater than 70 dB and detect the weak signal. In comparison to the sub-Nyquist 
scheme, Figure 19(b) shows a receiver using an 8 bit ADC with Nyquist sampling rate, the 
third input signal is too small to detect. 



Frequency (MHz) 



Frequency (MHz) 


Figure 19 . SFDR Test for Two Different Sampling Schemes 
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V. CONCLUSION 


This thesis has demonstrated how to utilize a noise source to reduce intermodulation 
distortion and increase the spurious free dynamic range of a channelized digital ESM receiver. 
This requires a wide bandwidth with resulting low dynamic range ADC to be used. In 
addition a sub-Nyquist sampling rate scheme was investigated which allows the full receiver 
bandwidth to be achieved with a smaller bandwidth but high dynamic range ADC. The data 
generated in this thesis are based on the selected frequencies of the input signals. Changing 
the frequency of the input signals results in slight variations of output data, but the principles 
remain the same. 

It is desirable to build a receiver with a wide bandwidth (greater than 1 GHz) and a 
high dynamic range (above 70 dB) in order to process multiple signals at the same time. 
Limited by the current ADC technology, a channelized digital ESM receiver meeting these 
specifications appears to be beyond the state-of-the-art. 

Although some spectrum estimation schemes such as sub-Nyquist rate sampling can 
achieve a wide bandwidth, it is unlikely to be used in real-time processing since they require 
intensive computation. It is difficult to achieve a wide bandwidth and high dynamic range in 
an ADC at the same time, but ADC technology is advancing at an astonishing speed, resulting 
in the digital receiver being a major trend of the future. 
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APPENDIX A 


% TITLE; NYQUIST SAMPLING RATE SIMULATION FOR ONE INPUT SIGNAL 
% SUBJECT: SIMULATION OF DIGITAL ESM RECEIVER SPECTRUM FOR ONE 
% INPUT SIGNAL USING FFT WITH TWO CHANNEL APPROACH 

% PROGRAM: FIG06.M 


clg; 

clear; 

fo=1.25e9; 

fs=0.5e9; 

fl=1.25e9; 

ifbw=0.5e9; 

nfft=512; 

N=8; 


% LOCAL OSCILLATOR FREQUENCY 
% SAMPLING FREQUENCY FOR IF SIGNAL 
% SAMPLE SIGNAL FREQUENCY #1 
% IF BANDWIDTH 
% NUMBER OF FFT FOR IF BAND 
% QUANTIZING BIT(S) 


% GENERATE I(t) and Q(t) SIGNALS 


dt=l/fs;. 

t=0:dt:(nfft-l)*dt; 
fl=round(fl *nfft/fs)*(fs/nfft); 
It=cos(2*pi*(fl-fo)*t); 
Qt=sin(2*pi*(fl-fo)*t); 


% SAMPLING TIME INTERVAL 
% SAMPLING PERIOD 
% MATCHING INPUT FREQUENCY 
% I CHANNEL SIGNALS 
% Q CHANNEL SIGNALS 


% QUANTIZATION PROCESS 

M=2^(N-1); % QUANTIZATION LEVELS 

for I=l:nfft; 
xx(I)=M*It(I); 
yy(I)=M*Qt(I); 
if xx(I)>0; 
x(I)=ceil(xx(I))/M; 
else x(I)=floor(xx(I))/M; 
end; 

if yy(I)>0; 
y(I)=ceil(yy(I))/M; 
else y(I)=floor(yy(I))/M; 
end; 
end, 

si^+j*y; % QUANTIZED TWO CHANNEL IF SIGNALS 

% FIND FOURIER SPECTRUM OF THE IF SIGNALS 
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SIF=fft(sif,nfft); % FFT OF QUANTIZED IF SIGNALS 

SIF=fftshift(SIF); % MOVE ZEROTH LAG TO CENTER SPECTRUM 

pif^SEF.*conj(SIF)/nfft; % PSD OF QUANTIZED IF SIGNALS 
PIF=pif7max(pif); % NORMALIZED IF PSD OF SIGNALS 

df^(fs/l e6)/nfft; % FREQUENCY SPACING OF FFT POINTS 

fx=0:df:(nfft-l)*df; % SET INPUT SPECTRUM BANDWIDTH 

% IF STAGE SPURIOUS SIGNALS 

pifii=sort(PIF); 

pifh(nfil)=[]; 

snr=10*logl0(sum(pifii)); % SIGNAL TO NOISE RATIO 

di-10*logl0(max(pifii)); % DYNAMIC RANGE 

subplot(211 ),plot(6c, 10*log 10(PIF),grid; 

text(0,26,'(a) IF Spectrum of 8 Bits ADC for One Sample Signal'); 

xlabel('Frequency (MHz)'),ylabel('Amplitude (dBc)'); 

axis([0 ifbw/le6 -100 20]); 

text(5 ,-3 0, ['SNR-, num2str(-snr),'dB']); 

% CALCULATIONS OF SFDR AND SINAD FOR DIFFERENT NUMBER OF BITS 

bits=12; 

forN=l:bits; 

M=2^(N-1); % QUANTIZATION LEVELS 

for I=l:nfft; 
xx(I)=M*It(I); 
yy(I)=M*Qt(I); 
if xx(I)>0; 
x(I)=ceiI(xx(I))/M; 
else x(I)=floor(xx(I))/M; 
end; 

ifyy(I)>0; 
y(I)=ceil(yy(I))/M; 
else y(I)=floor(yy(I))/M; 
end; 
end; 

sif=x+j*y; 

% FIND FOURIER SPECTRUM OF THE IF SIGNALS 

SIF=fft(sif); % FFT OF QUANTIZED IF SIGNALS 

SIF=ffishift(SIF); % MOVE ZEROTH LAG TO CENTER SPECTRUM 
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pif=SIF. *conj(SIF)/nffi; % PSD OF QUANTIZED IF SIGNALS 
PF=pif'max(pif); % NORMALIZED IF PSD OF SIGNALS 


% IF STAGE SPURIOUS SIGNALS 

ifs=sort(PIF); 

ifs(nffi)=[]; 

snr(N)=-10*log 10(sum(ifs)); 
snrt(N)=1.76+6.02*N; 
end; 


N=l:bits; 

subplot(212),plot(N,snrt,'o',N,snr,'*'),grid; 
xlabelCADC Bits'),ylabel('SNR (dB)'); 
text(6.2,12,'o; Theoretical Value'); 
text(6.2,5,'*: Experiment Value'); 

text(0,86,'(b) SNR for Different Number of Bits ADC with One Sample Signal'); 
fprintf('(a) SNR for Sample Signal =>\n'); 
for N=1 :bits; 

fyrintf('%g Bits ADC',N); 

4jrintf('SNR = %1.2f dB\n',snr(N)); 
end; 

fprintf('\n'); 
axis([0 bits 0 80]); 
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APPENDIX B 


% TITLE; NYQUIST SAMPLING RATE SIMULATION FOR TWO INPUT % 
SIGNALS 

% SUBJECT; SIMULATION OF DIGITAL ESM RECEIVER SPECTRUM FOR TWO 
% INPUT SIGNALS USING FFT WITH TWO CHANNEL APPROACH 

% PROGRAM; FIG07.M 

clg; 
clear; 
num=2; 
fo=L25e9; 
fs=0.5e9; 
fl=L23e9; 
f2=L27e9; 
ifb'w=0.5e9; 
nffi=512; 

N=8; . 

% GENERATE I(t) and Q(t) 

dt=l/fs; % SAMPLING TIME INTERVAL 

t=0;dt;(nfil-l)*dt; % SAMPLING PERIOD 

fl =round(fl *nffl/fs)*(fs/nffi); % MATCHING INPUT FREQUENCY 

f2=round(f2*nfft/fs)*(fs/nfft); % WITH FFT RESOLUTION 

It=cos(2*pi*(fl-fo)*t)+cos(2*pi*(f2-fo)*t); 

Qt=sin(2 *pi * (fl -fo) *t)+sin(2 *pi * (f2-fo) * t); 


% NUMBER OF SIGNALS 
% LOCAL OSCILLATOR FREQUENCY 
% SAMPLING FREQUENCY FOR IF SIGNAL 
% SAMPLE SIGNAL FREQUENCY #1 
% SAMPLE SIGNAL FREQUENCY #2 
% IF BANDWIDTH 
% NUMBER OF FFT FOR IF BAND 
% QUANTIZING BIT(S) 

SIGNALS 


% QUANTIZATION PROCESS 

M=2^(N-1); % QUANTIZATION LEVELS 

for I=l;nfft; 
xx(I)=(M*It(I))/num; 
yy(I)=(M*Qt(I))/num; 
if xx(I)>0; 
x(I)=ceil(xx(I))/M; 
else x(I)=fIoor(xx(I))/M; 
end; 

if yy(I)>0; 
y(I)=ceil(yy(I))/M; 
else y(I)=fIoor(yy(I))/M; 
end; 
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end; 

sif^+j *y; % QUANTIZED TWO CHANNEL IF SIGNALS 

% FIND FOURIER SPECTRUM OF THE IF SIGNALS 

Sff=fft(sif,nfft); % FFT OF QUANTIZED IF SIGNALS 

SIF=ffishift(SIF); % MOVE ZEROTH LAG TO CENTER SPECTRUM 

pif^SIF.*conj(SIF)/nffi; % PSD OF QUANTIZED IF SIGNALS 

PIF=pifi'max(pif); % NORMALIZED IF PSD OF SIGNALS 

d^(fs/l e6)/nfft; % FREQUENCY SPACING OF FFT POINTS 

fic=0:df;(nfft-1 )*df; % SET INPUT SPECTRUM BANDWIDTH 

% IF STAGE SPURIOUS SIGNALS 

pifii=sort(PIF); 

signal=pifh(nfft-num+I:nfft) % SIGNAL POWER 

pifh(nffi-num+l ;nfR)=[]; 

IMD=I 0*logl 0(sum(pifh(nfft-num-l O+I :nfft-num))) 

sinad=10*logl0(sum(signal)/sum(pifh));% SIGNAL TO NOISE AND DISTORTIONS 
sfdr=10*logl0(max(pifh)); % SPURIOUS FREE DYNAMIC RANGE 

subplot(21 l),plot(fx,10*logl0(PIF)),grid; 

text(0,26,'(a) IF Spectrum of 8 Bits ADC for Two Sample Signals'); 

xlabel('Frequency (MHz)'),ylabel('Amplitude (dBc)'); 

axis([0 fs/le6 -100 20]); 

text(405,-30,['SINAD=',num2str(sinad),’dB']); 

text(10,-30,['SFDR=',num2str(-sfdr),’dB']); 

text(310,10,'Input Signals','sc'); 

text( 110,10,'lntermods','sc'); 

% CALCULATIONS OF SFDR AND SINAD FOR DIFFERENT NUMBER OF BITS 

bits=12; 
for N=l;bits; 

M=2^(N-1); % QUANTIZATION LEVELS 

for I=l;nfft; 
xx(I)=(MHt(I))/num; 
yy(I)=(M*Qt(I))/num, 
ifxx(I)>0; 
x(I)=ceil(xx(I))/M; 
else x(I)=floor(xx(I))/M; 
end; 

if yy(I)>0; 
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y(I)=ceil(yy(I))/M; 
else y(I)=floor(yy(I))/M; 
end; 
end; 

sif=x+j*y; 

% FIND FOURIER SPECTRUM OF THE IF SIGNALS 

SIF=fft(sif); % FFT OF QUANTIZED IF SIGNALS 

SIF=fftshift(SIF); % MOVE ZEROTH LAG TO CENTER SPECTRUM 

pif^SIF.*conj(SIF)/nffi; % PSD OF QUANTIZED IF SIGNALS 
PIF=pifi'max(pif); % NORMALIZED IF PSD OF SIGNALS 

% IF STAGE SPURIOUS SIGNALS 

ifs=sort(PIF); 

signal=ifs(nffi-num+l:nfft); % SIGNAL POWER 

ifs(nffl-num+1: nflft)=[]; 
sinad(N)= 10 * log 10(suni(signal)/suni(ifs)); 
snrt(N)=I.76+6.02*N; 
end; 

N=l;bits; 

subplot(212),plot(N,snrt,'o',N,sinad,'*'),grid; 
xlabelCADC Bits’),ylabelCSINAD (dB)'); 
text(6.2,12,'o: Theoretical Value'); 
text(6.2,5,'*; Experiment Value'); 

text(0,86,'(b) SINAD for Different Number of Bits ADC with Two Sample Signals'); 

fprintf('(a) SINAD for Sample Signal =>\n'); 

forN=l:bits; 

fprintf('%g Bits ADC ',N); 
fprintf('SINAD = %1.2f dB\n’,sinad(N)); 
end; 

fprintf('\n'); 
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APPENDIX C 


% TITLE: SUB-NYQUIST SAMPLING RATE SIMULATION FOR TWO 
% INPUT SIGNALS 

% SUBJECT: SIMULATION OF DIGITAL ESM RECEIVER SPECTRUM 
% USING SUB-NYQUIST RATE WITH TWO CHANNEL APPROACH 

% PROGRAM: FIG18 M 


clear; 

num=2; 

fo=1.25e9; 

fsl=0.1e9; 

fs2=0.5e9; 

fl=L23e9; 

12=1.27e9; 

ratio=l; 

Tao=0.1e-9; 

nfft=512; 

N=12; 


% NUMBER OF SIGNALS 
% LOCAL OSCILLATOR FREQUENCY 
% SUB-NYQUIST SAMPLING FREQUENCY 
% NYQUIST SAMPLING FREQUENCY 
% SAMPLE SIGNAL FREQUENCY #1 
% SAMPLE SIGNAL FREQUENCY #2 
% RATIO OF SIGNAL #1 AND #2 
% DELAY TIME FOR DELAYED CHANNEL 
% NUMBER OF FFT 
% QUANTIZING BIT(S) 


% GENERATE EXTERNAL NOISE 


load noise.dat; 

Noise=0.001 *noise'; % SET NOISE LEVEL 

ifnffl<2048 
Noise(nffl+l:2048)=[]; 
end; 

% PART I: SIMULATION OF ADC USING SUB-NYQUIST SAMPLING RATE 

dt=l/fsl; % SAMPLING TIME INTERVAL 

t=dt:dt:nffi*dt; % SAMPLING PERIOD 

td=dt-Tao:dt:nfft*dt-Tao; % SAMPLING PERIOD FOR DELAYED CHANNEL 


s(t) 


% GENERATE INPUT SIGNALS 

fl=round(fl *nffi/fsl)*(fsl/nffi); 

t2=round(f2 *nfft/fs 1) * (fs 1 /nffi); 

stl=sin(2*pi*fl*t); 

st2=ratio*sin(2*pi*f2*t); 

st=stl+st2; 

std 1 =sin(2*pi*fl *td); 


% MATCHING INPUT FREQUENCY 
% WITH FFT RESOLUTION 
% SINUSOIDAL INPUT SIGNAL #1 
% SINUSOIDAL INPUT SIGNAL #2 
% SINUSOIDAL INPUT SIGNALS 
% SINUSOIDAL INPUT SIGNAL #1 
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std2=ratio*sin(2*pi*f2*td); % SINUSOIDAL INPUT SIGNAL #2 

std=std 1 +std2; % DELAYED SINUSOIDAL INPUT SIGNALS 

SIGNAL=sum(stl .^2)+sum(st2.'^2); 

NOISE=sum(Noise.^2); 

SNR=10*logl 0(SIGNAL/NOISE) 

It=cos(2*pi*(fl-fo)*t)+ratio*cos(2*pi*(f2-fo)*t); % UNDELAYED I CHANNEL 
Qt=sin(2*pi*(fl-fo)*t)+ratio*sin(2*pi*(f2-fo)*t); % UNDELAYED Q CHANNEL 
Itd=cos(2*pi*(fl-fo)*td)+ratio*cos(2*pi*(f2-fo)*td); % DELAYED I CHANNEL 
Qtd=sin(2*pi*(fl-fo)*td)+ratio*sin(2*pi*(f2-fo)*td); % DELAYED Q CHANNEL 

% QUANTIZATION PROCESS 


M=2^(N-1); % QUANTIZATION LEVELS 

for I=l:nffi; 
xx(I)=(M*It(I))/num; 
yy(I)=(M*Qt(I))/num; 
xxd(I)=(M*Itd(I))/num; 
yyd(I)=(M*Qtd(I))/num; 
if xx(I)>0; 
x(I)=ceil(xx(I))/M; 
else x(I)=^oor(xx(I))/M; 
end; 

if yy(I)>0; 
y(I)=ceil(yy(I))/M; 
else y(I)=floor(yy(I))/M; 
end; 

if xxd(I)>0; 
xd(I)=ceil(xxd(I))/M; 
else xd(I)=floor(xxd(I))/M; 
end; 

if yy(I)>0; 

yd(I)=ceil(yyd(I))/M; 
else yd(I)=floor(yyd(I))/M; 
end; 
end; 

sif=x+j*y; % QUANTIZED UNDELAYED CHANNEL IF SIGNALS 

sifd=xd+j *yd; % QUANTIZED DELAYED CHANNEL IF SIGNALS 

% FIND FOURIER SPECTRUM OF THE IF SIGNALS 

SIF=ffi(sif,nffi); % FFT OF QUANTIZED IF SIGNALS 

SIFd=ffi(sifd,nffi); % FFT OF QUANTIZED IF SIGNALS 

psd=abs(SIF); % OUTPUT AMPLITUDE 
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P SD=p sd/max(psd); 

df=(fsl/le6)/nfft; 

fx=0:df:(nfft-l)*df; 


% NORMALIZED IF PSD OF SIGNALS 
% FREQUENCY SPACING OF FFT POINTS 
% SET INPUT SPECTRUM BANDWIDTH 


subplot(211 ),plot(fx,20*logl 0(PSD)),grid; 

text(0,8,'(a) Output Spectrum of a 12 Bits ADC Sampling at fs=100MHz Sub-Nyquist 
Rate','sc'); 

xlabel('Frequency (MHz)'),ylabel('Amplitude (dBc)'); 
axis([0fsl/le6-120 20]); 

% IF STAGE SPURIOUS SIGNALS 

PSDN=sort(PSD); 

PSDN(nfft-num+l :nffl)=[]; 

SFDR=.20*log 10(max(PSDN)); % MAX SPUR LEVEL 

No=l; 

for m=l:nfft; 
if PSD(m)>max(PSDN) 

Tu=atan(imag(SIF(m))/real(SIF(m))); 

T d=atan(imag(SIF d(m))/real(SIFd(m))); 

Freq=round(abs((Td-Tu)/(2*pi*Tao))he6); 
fprintf('Input Signal %g',No); 
fprintfC = %g MHz\n',Freq+le3); 
f(No)=Freq+le3; 

No=No+l; 

end; 

end; 

text(82,5,['fl=',nuni2str(f(l)),'MHz'],'sc'); 

text(82,-6,['f2=',num2str(f(2)),'MHz'],'sc'); 

text(2,5,['SFDR=',num2str(SFDR),'dB'],'sc'); 

% PART II: SIMULATION OF ADC USING NYQUIST SAMPLING RATE 
N=8; 

dt=l/fs2; % SAMPLING TIME INTERVAL 

t=0:dt:(nfft-l)*dt; % SAMPLING PERIOD 

% GENERATE INPUT SIGNALS s(t) 

fl =round(fl *nffi/fs2)*(fs2/nfft); % MATCHING INPUT FREQUENCY 

f2=round(f2*nffl/fs2)*(fs2/nffi); % WITH FFT RESOLUTION 

It=cos(2*pi*(fl-fo)*t)+ratio*cos(2*pi*(£2-fo)*t)+Noise; % I CHANNEL SIGNALS 
Qt=sin(2*pi*(fl-fo)*t)+ratio*sin(2*pi*(f2-fo)*t)+Noise; % Q CHANNEL SIGNALS 


47 


% QUANTIZATION PROCESS 


M=2^(N-1); % QUANTIZATION LEVELS 

for I=l:nffi; 
xx(I)=(M*It(I))/num; 
yy(I)=(M* Qt(I))/num; 
if xx(I)>0; 
x(I)=ceil(xx(I))/M; 
else x(I)^oor(xx(I))/M; 
end; 

if yy(I)>0; 
y(I)=ceil(yy(I))/M; 
else y(I)=floor(yy(I))/M; 
end; 
end; 

siT=x+j*y; % QUANTIZED TWO CHANNEL IF SIGNALS 


% FIND FOURIER SPECTRUM OF THE IF SIGNALS 


SIF=fft(sif,nffi); 
SIF=fftshift(SIF); 
pif=SIF. *conj(SIF)/nffl; 
PIF=pif7max(pif); 
df^(fs2/le6)/nfft; 
fx=0: df (nfft-1) * df; 


% FFT OF QUANTIZED IF SIGNALS 
% MOVE ZEROTH LAG TO CENTER SPECTRUM 
% PSD OF QUANTIZED IF SIGNALS 
% NORMALIZED IF PSD OF SIGNALS 
% FREQUENCY SPACING OF FFT POINTS 
% SET INPUT SPECTRUM BANDWIDTH 


subplot(212),plot(fx, 10*logI0(PIF)),grid; 

text(0,8,'(b) Output Spectrum of an 8 Bits ADC Sampling at fs=500MHz Nyquist 
RateVsc'); 

xlabel('Frequency (MHz)'),ylabel('Amplitude (dBc)'); 
axis([0 fs2/le6 -120 20]); 


% IF STAGE SPURIOUS SIGNALS 


pifn=sort(PIF); 
pifn(nfft-num+l :nfft)=[]; 

SFDR=-10 * log 10(max(pifn)); 
fl =round(fl /1 e7) * 10; 
f2=round(f2/le7)*10; 

text(10,6,['SFDR=',num2str(SFDR),'dB’],'sc'); 
text(410,5, ['fl =',num2str(fl),'MHz'],'sc'); 
text(41 0 ,- 6 , ['f2=',num2str(f2),'MHz'],'sc'); 
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